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Abstract 

We have updated our analysis of the 9-year WMAP data using the collection of polarization maps 
looking for the presence of additional evidence for a finite ’cosmic ray foreground’ for the CMB. 
We have given special attention to high Galactic latitudes, where the recent BICEP2 findings were 
reported. The method of examining the correlation with the observed gamma ray flux proposed in 
our earlier papers and applied to the polarization data shows that the foreground related to cosmic 
rays is still observed even at high Galactic altitudes and conclusions about gravitational waves are 
not yet secure. Theory has it that there is important information about inflationary gravitational 
waves in the fine structure of the CMB polarization properties (polarization vector and angle) and 
it is necessary to examine further the conclusions that can be gained from studies of the CMB 
maps, in view of the disturbing foreground effects. 
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INTRODUCTION 


The ’structure’ of the Cosmic Microwave Background (CMB) is a fertile area for the 
testing of cosmological models, but the elimination of foreground effects due to processes in 
the Galaxy, particularly in the Halo, is not a trivial problem. The processes relate mainly to 
photons due to Cosmic Ray (Cl^ electrons and radiation from dust. The frequency used in 
the CMB maps (Bennett et.al. jd] hereafter referred to as the WMAP) is that for which the 
foreground contributions are a minimum. It is interesting to note that at this frequency 
the sum of the components gives almost the 2.7 K temperature spectral shape. 

In our recent publication ([l^), which conhrmed results reported by us about ten years 
ago (j^) with the one-year only WMAP data, Tegmark, de Oliviera-Costa & Hamilton js], 
we have shown that there is strong evidence that the procedure of eliminating the known 
and identified foreground contributions is not perfect. We have found some contributions 
from large Galactic structures which have not been taken into account yet. 

The method we used was based on the correlation of the explicit cosmic ray origin effects 
measured by the EGRET satellite Esposito et ah j^, i.e. the gamma ray flux, and the CMB 
temperature characteristics. 

In the present work we are still using the EGRET Gamma Ray measurement but com¬ 
paring the high energy photon flux with the polarization characteristics of the microwave 
radiation registered by the WMAP experiment. 

Some remarks are necessary about the sources of the data. Although the CMB ’tem¬ 
perature’ maps used earlier were cleaned of foreground, to the best degree available, by the 
authors, those for polarization - as published - have not been. There is thus an inevitable 
foreground contribution. Turning to the gamma ray map, the EGRET measurements have 
been superseded by the Fermi-LAT data. However, a comparison of the two sky maps shows 
that they are very similar. The features in the maps which are relevant to the present anal¬ 
ysis relate largely to the distribution of CR-target material, which is in the form of atomic 
and (lumpy) molecular hydrogen. It is the interaction of CR with the gas which produces 
the measured gamma rays. Our hypothesis regarding the CMB foreground is that some, 
at least, will be due to CR-heated dust and electron synchrotron radiation, both of which 
are only indirectly connected with the gas. Thus, and in order to preserve continuity with 
previous work, we continue to adopt the EGRET maps. 
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THE ANALYSIS 


The sky map of CMB temperature fluctuations has been analysed in many different ways. 
The most important and significant is the search for the cosmological structures observed 
against the background of random noise. On the other hand the observed sky map at any 
frequency of the electromagnetic spectrum has obvious structures related to the Galactic 
’shadow’ or to Galactic sources. These structures make a foreground on the GMB sky and 
we need to identify them and then extract them. For the temperature maps the process 
of the foreground reduction gives the set of ’’cleaned” maps, which contains to a great 
extent the real ’cosmological’ GMB structures, but the process of cleaning is not obvious 
and there could be different strategies leading to slightly different results. This is a source of 
systematics in the further analysis of the GMB data and in turn the source of uncertainties in 
the resulting cosmological parameters. For the polarization case the vector character of the 
data makes the cleaning procedure much more complicated. In addition, the small degree of 
polarization makes the flux under consideration much harder to determine experimentally. 
All this is probably the reason that the set of polarization data of WMAP available to the 
public do not contain the ’cleaned maps’ with any known cleaning procedure. 

The analysis of the WMAP polarization data shows a strong enhancement of the intensity 
of radiation close to the Galactic disk as would be expected for all foregrounds. 

The same structure is seen for the gamma ray flux in the EGRET map. We supposed 
that the EGRET map (with point sources excluded) is a view of the possible foreground 
on GMB maps caused by the cosmic rays and Galactic matter structures. The correlation 
between the EGRET gamma ray flux and the temperature fluctuations measured in this 
direction by WMAP was used in our previous papers {d, 10|. We have followed this way 
of constructing the sky map of the respective correlation coefficients for the polarization 
parameters of WMAP. 

The polarization is measured as a set of Stoke’s parameters {U and Q). For each point 
(pixel) in the sky the number of independent observations is also given. The statistics 
(exposure) of different parts of the sky is not uniform and there are patterns relevant to the 


Q.' 


detector structures and measurement technique as mentioned already in [Ti.This is a factor 
which contributes slightly to the reported average values of polarization. The polarization 
intensity, /, even if there is a completely random polarized light registered, is dehned as the 
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square root of + U^), and has non-zero average values. The square of / fluctuates as 
the sum of the squares of two randomly distributed quantities. Each value of U and Q is 
the average value of a number of measurements, so it should in principle obey a normal, 
Gaussian, distributions with width proportional to the square root of the measurement 
number. All Stokes parameter values (/, U and Q) measured by WMAP depend on how 
long each fragment of the sky was observed. 

This is very important when one is looking for the small structures of polarized CMB. In 
the present paper, however, we concentrate on big and very big structures, as were the ones 
we found in the temperature fluctuation analysis. 

To diminish the effect of random fluctuations we used the technique of the ’running av¬ 
erage’ thereby smoothing the small scale effects. The smoothing of the WMAP polarization 
parameter maps was done by calculating the value of each parameter for each point on the 
high resolution WMAP, averaging all the points within a cone of radius 10° weighted with 
the exponential weight: exp[—(angular distance/3°)^]. 

Such smoothed Q and U maps were used to derive the polarization intensity I and 
the polarization angle 6 = 1/2 arctan(17/Q). We dehne the angle with respect to the N-S 
direction: value 0 represents the N-S direction and 90° is for the W-E direction. 

The procedure of running average was used to smooth the WMAP data for all frequency 
bands. Further analysis has been performed for each band separately, but also we add the 
fluxes of Q, V and W bands to produce a less fluctuating picture. The results and conclusions 
do not change using the sum, so below we present some details only for this case. 

The input maps are given in Fig HI 

The polarization map shows the emission from the Galactic disk, which is not excluded 
from the data sets. There is also visible the Equatorial Plane which could be partially an 
effect of the statistics discussed above. The angle map does not contain these well dehned 
and understood structures, but in general the angle map could be misleading. At every point 
the value of 6^ = 1/2 atan(f//Q) is obtained whatever the values of Q and U are, however 
small is the total polarization flux itself. Due to the uncertainties of the measurement for 
small values of U and Q the polarization direction is highly undetermined (small fluctuations 
of Q around the value of 0 from a positive to negative value change abruptly the direction 
from N-S to W-E). The maps for the angle should therefore be used with great care. 

Maps for both these quantities were then compared with the EGRET intensity map 
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FIG. 1: The WMAP data smoothed (see text) sky maps for polarization (a) and polarization angle 
(b). The angle equal to 0 degrees represents the N-S direction while 90° is the W-E direction. The 


maps are obtained using the software of the HEALPix package 


ix ^ 


shown in Fig. |2j For each pixel the correlation with gamma ray intensity was calcnlated as 
a linear correlation coefficient for all point in the neighbonrhood of the pixel dehned as the 
directions in the sky within the cone of radius of 10°. 


. (E,W,) - (E,)(W^> ^ 

^ (TeCTw 

where is the EGRET gamma ray flux and Wj - the polarization or the polarization 
angle towards the point j. Parenthesis (...) in Eq.([T]) means averaging over all points in the 
cone centred at the point j; cte and cxw are respectively the squared variances of the EGRET 
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FIG. 2: The EGRET gamma ray intensity map. 

gamma ray flux and the WMAP polarization signals - within this cone. 

The correlation coefficient maps are shown in Fig. |3]for the polarization intensity (a) and 
the polarization angle (b). The hrst shows obvious correlation with the disk emission from 
the EGRET map. The disk component is limited to small b values. Above b = 30° it seems 
to vanish. The map in Fig|3](b) shows similar patterns limited again to 6 < 30°. 

In the sky far from the Galactic Plane there are some pattern of quite large sizes seen 
in Figs. [T] and [3l They could come from random fluctuations (’cosmological variance’) or 
they could be related to some processes determined by big structures in the Universe, most 
likely in the vicinity of our Galaxy, or some Galactic structures, or even quite local Galactic 
structures. 

To study this we have to look hrst at the size spectra of the patterns. 


Spherical harmonics 


To examine this structures we can use the spherical harmonic analysis. It is known 
that the distribution of the signals on a sphere can be expanded as a linear combination of 
spherical harmonics: 

OO i 


£=0 m=—i 

The effect of hltering (by the smoothing procedure applied to the raw data) was studied in 
Wibig & Wolfendale jJ. No additional power resulting from the ’running average’ we used 
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FIG. 3: The WMAP data smoothed (see text) sky maps for polarization (a) and polarization angle 
(b). The angle equal to 0 degrees represents the N-S direction while 90° the W-E direction. 

is observed in the spectrum for values of ^ below the respective averaging size. 

Spectra of the polarization and polarization angles measured by WMAP are shown in 
Fig. m We compared them with the randomized sky map spectra obtained by shuffling pixels 
on the original WMAP maps. Such maps do not contain any information about the Galaxy 
and the position with respect to any real structures, so they can be used as an estimation 
of the background for our foreground analysis. As we can see in Fig. 01 for ^ greater than 20 
the effect of the running average suppressed possible existing forms in the CMB sky. The 
curious lowering of the angle spectrum even below randomly shuffled pixel sky spectra for 
£ 20 is not understood: it is a specihc feature of the map in Fig. [T]d. 

The correlation coefficients maps suggest the existence of non-random low ^ patterns, 
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FIG. 4: Smoothed harmonic spectra (shown as Ylm ^ 1)) calculated for the polarization 

(a) and polarization angle (b) of the maps shown in Fig(T] (solid lines) compared with the spectra 
obtained for three examples of ’randomized’ skies (see text) spectra . 
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FIG. 5: Smoothed harmonic spectra calculated for the linear correlation coefficient of the EGRET 
intensity flux and the polarization (a) and the polarization angle (b) of WMAP as shown in 
Fig{T] represented by the solid lines compared with the spectra obtained for three examples of 
’randomized’ skies spectra. 









too. The power spectra of maps in Fig. [3] are shown in Fig. [5l We can see that the 
effect of the structures exists in the correlation maps. The unusual behaviour of the angle 
spectra around t = 10 — 20 has disappeared. Polarization and polarization angle correlation 
coefficient spectra look quite similar, only for £ of a few the large Galactic disk related 
spectral power of the correlations for polarization value is slightly higher than for angle. 

To see if the effect is only that related to the well visible Galactic disk, or is extended over 
high altitudes, we examined the spectra calculated for the parts of the sky with different 
cuts on Galactic altitude b. A search for more complicated big patterns is almost impossible, 
mainly because of the uncertainties of the measured values of the polarization light. Instead 
of plotting the whole spectra, we sum up the values shown in FigsJH and E] in the low £ 
interval: from 3 to 10. We dehned the respective ’spectral power density’ as 


10 


f3-io(^cut) — 




cut. 


E E I/. 


m 1 2 


X 


' + 1 ) , 


(3) 


i=S m=—i 

where G(6cut) is the solid angle of the analyzed sky covering the region of |6| > beat- 


RESULTS FOR HIGH ALTITUDES 

Results of our calculations of fa-io are shown in Figs. [6] and [71 

As we can see for the polarization value (Fig. [6^.) there are structures existing to Galactic 
altitudes of about 60° and these structures are not random patterns. In the case of the 
polarization angle (Fig. [6b) there is an enhancement of the power very close to the Galactic 
disk, and then the spectral power for low i goes close to the random sky spectral power, but 
is still about 20% higher. We would like to examine if this is by chance, or is a real effect. 

This question can be answered, if we check if these structures are related to the Galactic 
structures observed in the gamma ray intensity map of the EGRET experiment. The re¬ 
spective sums of power in the low £ parts of the linear correlation coefficient for polarization 
and polarization angle are shown in Figs. [6] and [71 

Gomparing the power contained in the correlation coefficient spectra three features are 
seen: 


- the low £ spectral power is observed both in the polarization intensity and in the 
polarization angle. 
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FIG. 6: ’Spectral power density’ calculated for the polarization (a) and the polarization angle (b), 
respectively. The results of WMAP, represented by the solid lines, are compared with spectra 
obtained for three examples of ’randomized’ skies spectra. 




FIG. 7: ’Spectral power density’ for the linear correlation coefficient of EGRET ffux and the 
polarization (a) and the polarization angle (b) shown as in Fig.[6l The results of WMAP represented 
by the solid lines compared with the spectra obtained for three examples of ’randomized’ skies 
spectra. 
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FIG. 8: Sum of low i (< 10) spherical harmonics for the correlation map of the polarization (a) 
and the polarization angle (b) . 

- in both cases the 6cut dependence is very similar 

- the existence of the non-random structnres is observed in the whole sky, np to 6 = 75°. 

This could stand for evidence that the polarization intensity and the angle are related to 
the gamma ray flux rather than to the Galactic structure itself. 

We plot in Fig. [H] the sum of 10 of the first £ harmonic spectral components for the linear 
correlation coefficient of the EGRET flux and the polarization value and the polarization 
angle. The maps consists of a number of blobs of size 20° — 30° which are distributed all 
over the high b regions. 
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CONCLUSIONS 


We have examined the CMB WMAPs in some detail from the standpoint of searching 
for Galactic foreground effects. They appear to be there, at all Galactic latitudes and not 
just near the Galactic Plane.This conclusion follows from a comparison of the polarization 
maps with those of Gosmic Gamma Rays. The presence of a signihcant foreground in the 
latitude range —50 to —65 degrees, is such as to render the interpretation of the BIGEP-02 
polarization data in terms of inflationary gravitational waves untenable at the present 
time (see, e.g., Adam et ah [l|). However, the fact that the foreground does diminish to 
some extent with increasing latitude, |6|, suggests that a more sophisticated analysis might 
be trustworthy. 
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